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Abstract 



We consider the coulomb gas model on the upper half plane with different boundary 
conditions, namely Drichlct, Neuman and mixed. We related this model to SLE(k,/?) 
theories. We derive a set of conditions connecting the total charge of the coulomb gas, 
the boundary charges, the parameters k and p. Also we study a free fermion theory in 
presence of a boundary and show with the same methods that it would lead to logarithmic 
boundary changing operators. 
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1 Introduction 

Conformal field theories have found many applications in classification of phase transitions 
and critical phenomena in two dimensions and other areas such as string theory. In particular 
the minimal models introduced in reveal many exact solutions to various two dimensional 
phase transitions like Ising model at critical point or three critical Ising model and so on. 
These models were first considered on the whole plane, but as many surface phenomena 
are very interesting to analyze, boundary conformal field theory was soon developed j2J. 
Essentially it was shown that conformal field theory in the half plane with proper boundary 
condition (BC) could be mapped to a whole-plane conformal field theory, the price you have 
to pay is to insert image fields in the other half plane. The idea was then applied to different 
problems such as turbulence 

On the other hand, recently a new method to investigate the so called geometrical phase 
transitions has been developed, which were previously described by conformal field theories. 
The new method, Stochastic Lowener Evolution (SLE) [I] is a probabilistic approach to study 
scaling behavior of geometrical models. For a review see jSJ H3 |S| . SLE's can be simply 
stated as conformally covariant processes, defined on the upper half plane, which describe 
the evolution of random domains, called SLE hulls. These random domains represent critical 
clusters. 

The idea of SLE was first developed by Schramm He showed that under assumption 
of conformal invariance, the scaling limit of loop-erased random walks is SLE2. (SLE's are 
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parameterized by a real number k, and are abbreviated as SLE/j) Also he claimed without 
proof that SLFjq is the scaling limit of critical percolation. The claim was proved later on by 
Smirnov [S]. He showed that in the scaling limit of percolation, conformal invariance exists 
and also using the new technic proved Cardy's formula jlUl lllj . 

As in SLE we are dealing with conformal mappings and geometrical phase transitions, 
one expects that there should be a relation between SLE and conformal field theory. The 
existence of a boundary is an essential point in SLE, so one would expect it would be related 
to BCFT's. An explicit relationship between the two was discovered by Bauer and Bernard 
^3 El EE21 EE] and afterwards by They coupled SLE K to boundary conformal 

field theory with specific central charge depending on the parameter k. These CFT's live 
on the complement of SLE hulls in the upper half plane In such situations, boundary states 
emerge on the boundary conformal field theory. The good point is that these states are 
zero modes of the SLE K evolution, that is they are conserved in mean. This means that all 
components of these states are local martingales of SLE re and hence one is able to compute 
crossing probabilities in purely algebraic terms. The key point of these result is existence of 
null vectors, just an in the case of minimal models. It turns out that the Verma module in 
the CFT part should contain level 2 null state, though higher level null vectors have been 
considered, too [2011^. The level two null state obeys the equation 

(2^2 + ^)^=0, (1) 

where L_f c , s are the Virasoro operatores. This equation naturally leads to second order 
differential equations on correlators containing the null state which are the same ones obtained 
from stochastic approach. 

More recently a generalization of these theories has emerged: the so called SLE(k, p) 
|28[ I29j . The new theory describes random growing interfaces in a planar domain which 
have markovian property and conformal invariancy. In fact SLE(k, p) is the minimal way to 
generalize the original SLE while keeping self-similarity. We will come back to the properties 
of such theories in the third section. 

Cardy [23| extended the correspondence between SLE K and CFT, to the case of SLE(k, p). 
In this case the equation in the CFT part is a little bit different from equation it has an 
extra term which comes from a current density J: 

(2L^ + ^L 2 _ 1 -J^ 1 L^\h) = 0. (2) 

If a current density is already present in the theory, then the equation would have a more 
physical meaning. The example considered in j2jj is a free bosonic field with the action 
(g/An) J(d(p) 2 d 2 z. The boundary conditions are piecewise constant Drichlet; that is on real 
axis we have 

4>(x) = 2-7T ajH(x — Xj) (3) 
j 

where H(x) is the Heavyside step function. 

The effect of a jump of strength a could be replaced with a boundary condition changing 
operator (j) a . The state produced by operation of this field is the one which satisfies the 
equation (J2J. The correlation functions associated with such field can be defined in the 
following way: first find a classical solution satisfying the boundary condition (BC). Calling 
the solution (f) c i, we write <f> = <j) c i + (j)' , where (j)' is assumed to vanish on the boundary. 
The partition function is then calculated to be Z = Z a Z'. Note that this happens, because 
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the theory is free. In boundary conformal field theory (BCFT) one can think of Z a as the 
correlation function {Ylj^aj}- 

Let's turn back to the specific problem. The classical solution of equation of motion 
satisfying the boundary condition (J3]is 

(t>d = i^2 a i ^ n (( z ~ x i)/(* ~ x j))- ( 4 ) 
Hence the partition function and the correlators are obtained to be 

*.-(n*,)-n(^p*- < 5) 

\ j I j<k v 7 



In a similar way the expectation value of energy-momentum density in presence of such 
BC is found, which can be thought of the correlation of energy-momentum tensor with the 
boundary fields (/> a 's. Conformal Ward identity comes afterward with which Cardy has found 
the action Lq, L_2 and L 2 _ 1 on the boundary fields. The fields are found to be primary with 
scaling dimension is h{ = gaf. Going back to the correlation functions ©, one finds that the 
two results are compatible if one assumes that the sum of all jump vanishes. 

To connect all these to SLE, one should look if it is possible to satisfy either equation Q 
or @- The first one is satisfied if we assume k = 4 and af = a* 2 = l/4g, so that hi = 1/4. 
The other one which is related to SLE(k, p) appears if we loose the second condition, which 
leads to 

2L_ 2 <j) ai = 2L 2 _ l( j> ai - '— 3— L-i(f> ai , (6) 

3 3 1 

where pj = (aj — oti)(l — (ctj/a*) 2 ) and prime means that j is excluded from the summation. 
This equation can be written as @ in terms of J = —2ig 1 / 2 d<j) which is conserved due to 
equations of motion. 

Cardy then mentions that it is possible to take Neuman boundary condition (NBC) instead 
of Drichlet boundary condition (DBC), one just have to consider the dual description of the 
problem with NBC on dual fields <fi with insertions of vertex operators e to "^. 

In the next section, we will consider the coulomb gas action with a charge at infinity and 
investigate different boundary conditions. The next chapter is devoted to the connections 
to SLE and SLE(k,/j). In the last section we will apply the same method to the fermionic 
action of c = —2 and find logarithmic operators and some structures which has been found 
before [251 • 

2 Boundary Operators in Coulomb Gas Model 

Coulomb gas is defined via the following action: 

S = ^j[g{dv? + 2QR V ] (7) 

where R is the scaler curvature associated with the background metric. We are interested 
in the case where <p is defined on the upper half plane. The boundary conditions can have 
several forms: Dirichlet (DBC), Neuman (NBC) or mixed boundary condition, but not all 
of them are conformally invariant. A boundary condition is conformally invariant if the 
antiholomorphic part of the energy momentum tensor in the lower half plane is the analytic 
continuation of its holomorphic part in the upper half plane. This is satisfied if T xy = at 
y = where 

T xy = ^d x <pd y (p - ®d x d y cp. (8) 



3 



In NBC case, where dy(p\ y= o = 0, the above condition is automatically satisfied, while in DBC 
case, it is not satisfied unless if Q = 0. Let's take Neuman boundary condition, we construct 
our desired boundary condition in the following way: suppose we have some electrical charges 
at points Xj located at the boundary. We would like to find the boundary changing operators 
associated with these electric charges. The classical solution which satisfies Laplace equation 
together with the above property has the form 

In this solution the normal derivative of field vanishes and hence the solution satisfies NBC. 
For consistency, (p has to behave like —2 Q ln(zz) when z, z — > oo, so one should assume the 
condition J2 Xj = —2Qg. 

the next step is to calculate the partition function. If we write ip = tp c \ + cp', with the 
condition d y ip'\ y= o = 0, we'll have S[(p] = S[(p c i] + S[np'], and the partition function associated 
with S[ipci] is derived to be 

j<k v a / 

Here a is the UV cut off. Let Z' be the partition function associated with iSfy?'], then we 
define the boundary correlation function to be 

<n<M* J )> = § (ii) 

3 

Also any expectation value (O) in presence of this boundary condition can be defined in 
the following way: 

{0) ^ (12) 

In particular the expectation value of energy-momentum tensor T can be computed to give: 

(T) = - 9 (dv d ? + 2Qd\ d = \ £ ^ + 2^ ^ 

9 ~£ \ z ~ x i)\ z ~ x k) 9j\z-*jY 

Taking the limit of z — >• Xi one can obtain the conformal weights of the operators <p(xi), 
h\i = Aj(Aj + 2Q)/g. The result suggests that these operators are related to vertex operators 
exp(iA^) in the coulomb gas model. 

Now let's investigate the Dirichlet boundary condition with Q = 0. The field is taken 
to be piecewise constant with some jumps 2na at the points Xj's. In this case, the classical 
solution of the equation of motion is 

^=^Y,^H Z —^) (14) 

Z Xj 

Note that in the presence of external charge, second term of offdiagonal energy-momentum 
tensor (|SJ) does not vanish at the points Xj. 

At this level, there are no conditions on ay's, and Z c i can be obtained as in the case of 
NBC. The result is just the same as equation (|1U|). you just need to transform g — ► 1/g and 
Xj — > ctj. Similar to the NBC case, there appear some operators on the points Xj. Also the 
definition of correlation functions of such fields is the same as equation It can be shown 
that these operators are related to magnetic vertices. 
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Again one can define the expectation values of any operator via a relation similar to 
equation (|12|). The related fields are then found to be primary ones with conformal weight 
h a . = get?. Now if we look back to the correlation functions of f a ^s, one observes that they 
are consistent with the weight derived, if the sum of all jumps add up to zero, that is the 
value of cp on both far ends of the boundary should be the same[2*7|. 

The other BC is the mixed one (with Q = 0). In this case (p c [ has the following form: 

V<* = J2~J ln (z ~ ~ Xj) + iJ2 a i ln f f _ V y ^j ( 15 ) 

The above BC means that some electric and magnetic charges are fixed on points Xj and 
Uj. It turns out that in the classical action, the two parts related to magnetic and electric 
charges decouple. The results are a combination of the previous two cases. If Xj = yj there 
is a electric-magnetic vertex operator O e ^ m with the scaling dimension h = (^/got + X/^g) 2 . 
It can be argued that for consistency we must have ^j a j = in addition to previous 
conditions. 

As in SLE, one can use similar bulk operators, rather than the boundary ones, in brief 
we will talk about bulk correlation functions. In the case of DBC one is able to use methods 
similar to the image method in ordinary electromagnetic theory, so if there is a magnetic 
charge at Zj then another magnetic charge with opposite sign should be located at Zj. Thus 
the neutrality condition is satisfied automatically. But in Neuman case, the image charge 
has the same sign and then the neutrality condition is not satisfied automatically. The two 
point functions of Dirichlet and Neuman BC's have the following form 



\z — wr\z — w* 



(O mi O m2 ) — \m 2 (j \ 

(Imz) i (Imw) 



z — w 



z — 



(17) 



These correlation functions could be used to develop a SLE(k, p) theory. 



3 Relation to SLE(k, p) 

In this section we would like to relate the result derived in the previous section to the new 
and interesting topic SLE(k, p). As we mentioned in introduction, SLE(k, p) is somehow a 
more general case of SLE K which describes random growing interfaces in a planar domain 
which have markovian property and conformal invariancy. This evolution has the following 
form 

2 n 

dt9t(z) = , g (z)=z, dW t = V^B t + y £ w ( V (18) 

g t {z)-W t iW t -gt(xj) 

where Xj's are some arbitrary points not necessarily on the boundary, and gt(z) is a conformal 
map from a subset of upper half plane, Ht, to the whole upper half plane, H. In other words, 
Ht = H\K — * H where Kt is the complement of Ht in the upper half plane. The subset Kt 
which is called the hull of evolution, is just the growing domains mentioned above. Setting 
all of the parameters pj to zero, one arrives at the ordinary SLE K , but with nonzero pj we 
may have several random curves growing from the points on the boundary [201 • SLE(k, p) 
has many interesting properties as you can see in |281 129j . 
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As we said earlier, Cardy [27| has established a connection between SLE(k, p) and a CFT 
with Drichlet boundary changing operators inserted at points Xj. He derives the properties 
of the boundary operators <pj which leads to equation © which is a specific case of the more 
general equation 

needed to establish a connection to SLE(k, p). 

Let's see how the different choices for boundary conditions in the previous section affect 
the derived SLE(/c, p). Taking g = 1, for the case of NBC, calculation of 2L_2 — f-^-i yields 



2Ln6, - 
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E ^^' ^p-H'lE ,,. -y fe (20) 

If we impose the two conditions A 2 = — = A 2 and Q = — - — j~ — ~i then 2L_2 — %I?L\ will be 
zero. Also we will have 

(6 - «)(3« - 8) 



^ = ^ (21) 



2k 1 2k 

which is similar to the case of SLE K . 

On the other hand, if we impose only the condition Q = — — — - we will have 

2L_ 2 </> - ^L 2 ^ = (2 - 2kA 2 ) £ -^-L^ (22 ) 

j,fe ^ 

which is the equation in SLE(k, p). 

There is an equivalent to study this problem. Bauer, Bernard and Kytola (BBK) have 
shown that SLE(k, p) could be obtained starting from correlation functions of boundary 
fields |30| . Later Kytola has applied this method to the coulomb gas model, that is he has 
considered the boundary correlation functions of a coulomb gas and derived the corresponding 
SLE(k, p) |31j. This method, although is a little bit different from the cardy's one, has the 
same consequences. By means of our results in the previous section, and using the method of 
BBK, one is able to derive the corresponding SLE(k, p). Namely the correlations (|10|) and the 
one obtained in other BC's are very suitable to do this. Additionally, Schramm and Wilson 
|32j and also Friedrich and Bauer [33]. argued that it not necessary to take the operators just 
on the boundary and one is able to construct SLE(k, Pbound.-, Pbulk) using both the operators 
on the boundary and in the bulk. To do this, the correlation functions ()16|) are helpful. The 
result are very similar to the ones we have derived, you should just replace Xj with Zj where 
Zj's are points in the upper half plane. 

Another point to mention is that coulomb gas action is not the only one to produce these 
results. Zamolodchikov and Zamolodchikov [31] have considered the Liouville action 

Sl= [ (^(d^f + pe 2 ^) (23) 



,47T 

with the boundary conditions 

4>(z,z) = -21ogzz + 0(l) \z\^oo (24) 

4>(z,z) = -2r]i\og(z - Xi)(z - Xi) + 0(1) \z\ -> Xi (25) 
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and have calculated the classical action to be 

Sl-cI = 2r liVj lo g \ x i ~ x j\ 2 ( 26 ) 
i<j 

where the condition rji = 1 should be imposed. As you see this is very similar to the 
correlation functions in the case of coulomb gas and will lead to the same SLE(k, p). 

Also Fateev, Zamolodchikov and Zamolodchikov [HE] have considered Liouville action 
on the disc with proper boundary conditions and derived the correlation functions of the 
operators in the bulk and on the boundary. Hence it should be related to a radial SLE(/t, p). 



4 Fermionic Boundary Changing Operators and LCFT's 

In this section we will apply the same method not to the action of a bosonic field, but to the 
action of a free fermionic fields: 

S = f 8689 (27) 



where 9 and 9 are grassman variables. This action is related to c = —2 conformal field 
theory which is known to be a logarithmic conformal field theory, that is, there can be 
found logarithmic terms in correlation functions of the fields inside the theory [22] • In such 
theories there exists pair of primary fields which transform into each other under conformal 
transformations. 

Before considering the action ()27|) we would like to review a powerful method developed 
to investigate LCFT's, namely the nilpotent weight method. The method is primarily based 
on a composite field, containing the two primary fields and a nilpotent variable which acts 
as a part of filed's weight f.24 a , that is, from the primary fields <j) and ip we construct the 
field <£(z,/3) = (ft(z) + (3ip(z) where (3 is a nilpotent variable. This fields acts as a primary 
field under conformal transformations with the weight h + 0, assuming (f) has weight h. The 
idea was then generalized in [2Sj, where the nilpotent variable was taken to be product of a 
grassman variable and its conjugate, hence the composite field contained four different fields, 
two bosonic and two fermionic. This structure was observed before in c = —2 theory by 
Kausch [2Hj- In this language the composite is written to be 

rj) = <j)(z) + fj£{z) + £(z)tj + fjrjip(z) (28) 

The fields 4> and ip are the same as previous ones and £ and £ are the new primary fermionic 
fields. The two point correlation function of such fields can be found exploiting conformal 
invariance: 

/*/ w x\ (m + m)(vi + m) / on s 
{H*i,mM*2,m)) = {zi _ Z2) 2 h+fjim+m2 ( 29 ) 

This equation could be read through expanding both sides of it, in terms of grassman variables 
to find correlation function of individual fields. 

With this brief review, we will go back to our problem, a c = —2 theory with different 
boundary conditions. Assume that the boundary condition is Drichlet with discontinuities otj 
at points Xj. We should do the same procedure as in the bosonic case. The classical solution 
satisfying this BC is 

^ = *E«,m(f^) (30) 

Here ay's are grassman variables in contrast with the bosonic case. The next step is to 
calculate the classical part of the partition function. As the solution ()3U|) is just the same 
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as the bosonic one, the result is the same, too. The only difference is that we have 6 in the 
action and this variable's jump is a rather than a. So we have 

*.=n(^r*- m 



a 



Now we can assign this partition function to correlation of boundary fields. Before doing this 
we have to impose a condition similar to charge neutrality in the bosonic case, that is we 
should assume that the sum of all the jumps has to be zero. To assure that this happens we 
define 

<n<>(«i>*j)> = z * x <HE«;)<HE«;)- ( 32 ) 

Here tfy's are the boundary fields and could be expanded in terms of ay's to have a) = 
(p(x) + a£(x) + £(x)a + aaip(z). Let's investigate the two point correlation of such fields. As 
the delta function of a grassman variable is itself, we'll have 

(i?(zi, ai) , &(x2, a 2 )) oc (x\ - X2 ^ a ^ +a ^) + a 2 )(ai + a 2 ) (33) 

Comparing with equation ()29|) and noting that sum of ay's is vanishing, we observe that we 
have found the correct result. Note that the weight of this field is zero as expected in c = —2 
theory. It may be possible to construct a non zero weight boundary field if we let the action 
(|27jl have a bosinic part. Higher correlation functions could be also derived in the same way, 
but now we would like to find the OPE of the fields with energy-momentum tensor. The 
expectation value of energy momentum tensor could be read easily through 



{T(z)) a = -Be d de d = 7 _ y < 34 ) 

jk *- x i Xk ' 



which leads to the OPE: 



= + — E' + + • ■ ■ . (35) 

which means that i? fields have weight equal to act, which is consistent with previous results 
|241 125 j . Using these boundary operators and their correlation functions, one is able to find 
the related SLE(k, p), exploiting the method introduced by Bernard et al. [3*U] . 
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